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Abstract 

We perform the stabihty analysis of the Schwarzschild black hole in f{R) gravity 
with the parity-violating Chern-Simons (CS) term coupled to a dynamical scalar field 
9. For this purpose, we transform the f{R) gravity into the scalar-tensor theory by 
introducing a scalaron (/>, providing the dynamical Chern-Simons modified gravity 
with two scalars. The perturbation equation for the scalar 9 is coupled to the odd- 
parity metric perturbation equation, providing a system of two coupled second order 
equations, while the scalaron is coupled to the even-parity perturbation equation. 
This implies that the CS coupling affects the Regge- Wheeler equation, while f{R) 
gravity does not affect the Zerilli equation. It turns out that the Schwarzschild black 
hole is stable against the external perturbations if the scalaron is free from the tachyon. 
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1 Introduction 



f{R) gravities [H El [3] have much attention as one of strong candidates for explaining the 
current accelerating universe [1]. f{R) gravities can be considered as Einstein gravity with 
an additional scalar (scalaron). For example, it was shown that the metric-/(i?) gravity is 
equivalent to the cjbd = Brans-Dicke theory with a certain potential [S]. 

On the other hand, the Chern-Simons (CS) modified gravity was obtained by adding 
a parity-violating CS term to the Einstein-Hilbert action, where the CS term couples to 
gravity via a CS scalar field 6 [B]. Originally the coupled scalar field 6 was considered 
as a prescribed function, but on later this choice was not regarded as the well-motivated 
one. Indeed, the dynamical Chern-Simons (DCS) modified gravity has been formulated by 
treating the scalar field ^ as a dynamical field [7j. For a review on the CS modified gravity, 
its astrophysical consequences, see [5] and for its critical gravity on the AdS4 spaceimes, 
see [H1[IU] 

It is very interesting to investigate the Schwarzschild black hole obtained from a modified 
gravity of the f(R) gravity with the dynamical CS term because astrophysical black holes 
are the most promising objects to probe the strongly gravitational field region of a modified 
gravity. The first study of the /(-R)-black hole stability has very recently been performed in 
the f{R, G) gravity [11]. In its scalar-tensor theory [12], the even-parity perturbations were 
affected by the scalaron and thus the black hole was stable against the whole perturbations if 
the scalaron did not have a tachyonic mass. In the context of the DCS modified gravity, the 
black hole perturbation has been carried out in [13], which indicates that if the background 
CS scalar ^ is a non-trivial, there was a serious mixing between odd-and even-parity metric 
perturbations. On the other hand, if ^ = or const., odd-and even-perturbations were 
decoupled as in Einstein gravity and odd-perturbations are affected only by the CS scalar 
field [13]. The odd-parity and CS scalar perturbations were described by a coupled system 
of two second-order equations, which has shown that the black hole is stable in the DCS 
modified gravity [T5] . 

Very recently, there was a perturbation study on the black hole in the context of /(-R, C) 
modified gravity with C the CS term [TB]. The black hole is unstable because the per- 
turbed Hamiltonian is not bounded from below, due to the CS term. In order to avoid the 
instability, either R = const, or = is required. In this case, number of physically 

propagating degrees of freedom are three, one from odd-parity and two from even-party and 
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scalaron because the f{R,C) modified gravity belongs to the non-dynamical CS modified 
gravity. Those modes are too strongly coupled to decouple three independent modes, which 
shows a distinctive feature of a parity-violating theory. However, the no-ghost condition of 
^^^bb"^ ^ no-tachyon condition of ^ 9^2''"^ > survive as in f{R) gravities. 

In this work, we wish to perform the stability analysis of the Schwarzschild black hole 
in f{R) gravity with the parity- violating CS term coupled to a dynamical scalar field 9. In 
order to avoid the difficulty with fourth-order derivative terms, we first transform the f{R) 
gravity into the scalar-tensor theory by introducing a scalaron 0. This will provide the DCS 
modified gravity with two scalars, which means that four modes are physically propagating 
degrees of freedom. Interestingly, the perturbation equation for the CS scalar 6 is coupled 
to the odd-parity metric perturbation equation, providing a system of two coupled second- 
order equations, while the scalaron is coupled to the even-parity perturbation equation. 
This enables us to perform the stability analysis of the Schwarzschild black hole obtained 
from /(i?)+DCS modified gravity theory completely. To make all things clear, we mention 
our notations. The metric signature is (— , -|-, +, +). The Riemann, Ricci tensor and Levi- 
Civita tensor are defined by 

1 



^^J.'- ucr 



A 



R 



R". 
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2 f{R) gravity with the DCS term 

Let us consider f{R) gravity with the dynamical Chern-Simons term in four dimensions 
which is given by 



1 

2^ 



f{R) + -*RR-aV,eV,e 



(2.1) 



where = SttG, a is a dimensional constant, and *RR = * R^^^'^ R^^^^^ is the Pontryagin 
density with 

1 



.p.- (2-2) 

Here e^'^'*^ denotes the four-dimensional Levi-Civita tensor. It is well known that the action 
can be rewritten by introducing a scalaron field as follows [T7] : 



1 
2^2 



(pR - 1/(0) + ^*RR - av^^ev^e 



(2.3) 
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with the potential V{(j)) = <f)A{(f)) — f{A{(f))). Note that the mass dimensions of 0, 6, and 
a are given by [(f)] = 0, [9] = —2, [a] = 4, respectively. Varying for the fields g^j,^, 0, and 9 
lead to the following equations: 

= -C,, + a (V,9VJ - ^ Vp^Ve) , (2.4) 

R = V'{4>), (2.5) 

V^9 = — -*RR (2.6) 

8a 



where ' denotes differentiation with respect to 0, and (7^,^ takes the form 

1 



C,. = V, 9 e''']y\,\R,y + - V,V. e e,'^'R^ (2.7) 



We take the trace of (12. 4p to rewrite (12.51) as the scalaron equation 

3V^0 + 2V{(f)) - (f)V'{(j)) = -2aV^9V^9. (2.8) 

Also we can express Eq. (l2.4p to be 

(t>R^.u - \g,.uV{(p) - ^g^uV^<p - V^Vu(t) = -C^. + aVf,9VJ. (2.9) 

Taking the restricted background valueil] as 

9 = const., = = const., 1/(0) = V'{^) = 0, ^"(0) ^ 0, (2.10) 

the solution to the Eqs. (12.41) . (12.51) and (12. 6|) is given by the Schwarzschild spacetime 

^^Sch ~~ g ^ydoc^ 

= -f{r)de + j^^+r\d^l + sm^^,d^l) (2.11) 



^ When taking these values, it gives the background spacetimes with the constant curvature scalar R 
which provides an easy step to find the solution of /(i?) gravity. In obtaining the constant curvature-black 
hole solutions (for example, Schwarzschild and Schwarzschild- (A) dS black holes), it seems that there is no 
difference between = and 9 = const.. Hence we choose 9 = const, here. Note that cf) corresponds to f'{R) 
in the original f{R) gravity and R = = 2y(0)/0 from Eqs. (|2.5p and (|2.8p . In this work, we focus on 

the Schwarzschild black hole solution with ^ = 0, which implies that V {(j)) = V{(}>) = 0. We mention that 
this is possible to occur when choosing a limited form of /(i?) gravity: /(i?) = aiR + a2R^ + - ■ ■ [H [H HO] ■ 
In this case, one finds that (f> = ai. 
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with the metric function 

/(r) = l-— . (2.12) 

r 

Now we introduce the perturbation around the background metric as 

Qtiu = (Jfiu + h^i,. (2.13) 
The perturbations around the background solution (p and 9 are given by 

6 = 6 + 66, (f) = ^ + 6(j). (2.14) 
The hnearized equation to (12.91) can be written by 

^6R,,ih) - ^V0y"(0)50 - V,VM = -SC^, (2.15) 
where the hnearized quantities of 5R^y{h), 6R{h), and 5C^y{h) take the forms 

6R{h) = V^V^h^^-V^h 

SC,Ah) = ^V,V. 5^ e(/^^i?%)^,. (2.16) 

In these expressions, the "overbar" denotes the background quantities. From Eq. fl2.8p and 
(12. 6p . we obtain the hnearized- scalar on equation 

V2 - l0V"(0)l50 = (2.17) 
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and the linearized-6' equation 



^'S6 = -^e^-^P'^R\^,V,V,hi (2.18) 

3 Perturbation analysis 

The metric perturbations h^^ are classified according to the transformation properties under 
parity, namely odd sector [ho, hi) and even sector {Hq, Hi, H2, K). However it is nontrivial 
task to show how the decoupling process goes with two scalar fields [56, 5(f)) wellj^ In order 



^It turns out that for 9 ^ const., there was mixing between odd and even modes in Chern-Simons 
modified gravity |13j . Here we can avoid this difficulty by choosing 9 = const.. 
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to see this explicitly, we must consider the full metric perturbation as 



''/if 



( Ho{r)Y H^{r)Y -|^/^o(r) sm^,d^Jh,{r) \ 

-S^^iW sin(^i(9<^,F/ii(r) 
r^YK{r) 

\ sin (/?i9(p^F/io(r) smipid^p^Yhiir) r"^ siv? ipiY K{r) J 

with Y = Y^^'^^{ifi, ip2) spherical harmonics. The form of 59 and 50 are given by 



Ho{r)Y 
Hi{r)Y 



H2{r)Y 



ikt 



(3.1) 



59 



Ye 



-ikt 



Mr] 



Ye 



-ikt 



(3.2) 



Substituting Eqs. fl3.ip and 03.21] into Eq. fl2.15p and after tedious manipulations, we find 
the perturbation equations for ten components as 



(t,t); 

(r,r); 



^-ikt^^Y = 



e-'^^E2Y = 



-ikt 



-ikt 



-ikt 



-ikt 



-ikt 



-ikt 



-ikt 



-ikt 



Esd^.Y + Oid^,Y 

E^d^,Y + O^d^.Y 
E^Y = 
E,d^,Y + O^d^^Y 

E,d^,Y + 0,d^,Y 











EqY + Ejd^J + 0,d^,Y + o,d^,d^,Y 



Esd^,Y + Ejd^.d^^Y + OjY + Osd'^Jj = 
EgY + Erdlj + Eiod^,Y + Ogd^.Y + Oiod^,d^,Y) = 0, (3.3) 
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where Ei with i = 1, ■ ■ ■ ,10 are functions of {Hq, ifi, i72, K, $) and Oi with i = 1, ■ 
are functions of (/iq, hi, ip) (see Appendix for the details). It is important to note that for 
L > 1, the perturbation equations fl3.3p imply twenty conditions like 



Ei = 0, Oi = 0, for z = 1, ■ ■ • , 10 



(3.4) 



which mean that ten perturbation equations can be decoupled into two classes: odd-parity 
{{Oi}) and even-parity {{Ei}). 
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For the even-parity case, we observe that the condition of = yields 

Ho{r) - fH^ir) - §^<l>(r) = 0. 

(pr 



(3.5) 



By using the above condition together with Ei = {i = 1, ■ ■ ■ ,6), one finds the central 
constraint equation as 

{Xr' - 2 + rrV'} Ho + |2A;V/-i + 2/ + rf + ^-^r\f'f - K 



X 



2ik 



2tkr + - 2A - 4/ - 2k'r'f-' - -r\fY — = 



0r 



(3.6) 



where A = L{L + 1). Manipulating two equations of -E2 = and -Es = by using the 
Eqs.dSS]) and ([SSD lead to 



d_ 

dr 



K 



0r 



2r(A-2/ + r/') 



K 



2A/-4A;V-A' (Hi 



r J 2zr2(A - 2/ + r/') V k 



^ V"^; " i{Xf-2P + rff') V ^ + ^ 

3A/-7V2 - 2/' + r/-7" - 2A;V/-i 

k 



X-2f + rf' 

Now we introduce the tortoise coordinate (r* = / y) and a new field defined by 



(3.7) 



(3.^ 



M 



1 



pq — h 



0r 



k 



(3.9) 



where 



g(r) 



p(r) 



A(A + ly + 3AMr + 6M^ , , , i(-Xr^ + 3AMr + 3M2) 

— , h{r) = 



r2(Ar + 3M) 
„2 ~ A 
r-2M' ^^2"^' 



(r-2M)(Ar + 3M) 



As a result, from the Eqs. (13.71) . (13.81) and (13. 9p we arrive at the Zerilli equation 

d^M 



dr*"^ 



+ 



k'-Vz 



where the Zerilli potential is given by 

2A2(A + l)r3 + GA^Mr^ + 18AMV + ISM^ 



Vz{r) = f 



r3(Ar + 3M)2 



(3.10) 



(3.11) 



(3.12) 
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The potential Vz{r*) is always positive for whole range of — oo < t* < oo , which implies 
that the even-parity perturbation is stable, even though the scalaron $ is coupled to making 
the even-perturbation [12]. In addition, using the tortoise coordinate (r*), the scalaron 
equation (12.171) becomes 



dr 



*2 



0, 



where the scalaron potential V$ is given by 

A 2M 



/(- 



+ 



(3.13) 



(3.14) 



with m| = (f)V " {(f)) / "i . The potential V$ is always positive exterior the event horizon if the 



mass squared m| is positive |12]J^I 

On the other hand, for odd-parity perturbation {{Oi = 0}), the first five equations 
provide three: 



Oi = or O2 = 



r\-4:M + Xr)ho - rf(2kir^ - ^M^p + kir^h\ + ^Mr^/;' + r^h'^ = 0,(3.15) 

\ (f) (f) / 

O3 = or O4 = 

(3.16) 
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-ikr^ (2ho - ikrhi - rh'^ + rV(A - 2)hi + -^ikM^jj = 0, 



O, = 



ikr^ho - (2M - r)|2M/ii - (2M - r)rh[^ = 



(3.17) 



and all remaining equations Oi with i = 6, - ■ ■ ,10 are redundant. Introducing the tortoise 
coordinate and a new field Q defined hj Q = fhi/r, the above three equations (13.151) ^ (13.171) 
become one coupled second-order equation 



d^ 



dr 



*2 



Q+<k'-f 



A 6M\ 



) Q 



QikMf 



(3.18) 



where Q = (j)Q. Also, the perturbation equation (I2.18P for the dynamical scalar 9 becomes 
a coupled second-order equation 



d'' 



dr 



*2 



18Af^ 

r^a 



■SM^^^e. (3.19, 



3ln the original f{R) gravity, the quantity of <?!)V^"((/))/3 corresponds to /'(0)/3/"(0) [0 4^ /'(O), F" ^ 
l//"(0)]. Therefore, the condition of m| > imphes no-tachyon (/"(O) > 0) if /'(O) > (no-ghost) in 
f{R) gravity. 
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where a = (pa. This is an important feature of CS couphng to f{R) gravity. Actually these 
coupled equation are the same found in [TB] . Hence, it is clear that the black hole is stable 
against the perturbations of Q and ip when using two independent numerical approaches 
of time evolution and a formation of frequency domain employed in Ref.|15]. 

Finally, we wish to mention the /(i?)-form dependence on the stability of the Schwarzschild 
black hole. In writing down two Eqs. fl3.18p and fl3.19p . we introduce two new variables 
Q = (f)Q and a = (pa which show the connection to the original /(-R) gravity because 
= /'(O). As was mentioned in footnote 1, our analysis is valid for a limited form of 
f{R) = aiR + a2R^ + ■ ■ ■ . In this limit from, we have /'(O) = ai = 0, which is fixed by 
choosing the limited f{R) gravity. In general, we can say that different f{R) theories with 
different corresponding (p have different model parameters. However, as far as the constant 
curvature-black hole stability is concerned, we expect that Eqs.f l3.18|) and fl3.19p remain 
unchanged except (p, leading to the stable Schwarzschild black hole. 

4 Discussions 

In this work, we have performed the stability analysis of the Schwarzschild black hole in 
f{R) gravity with the parity-violating CS term coupled to a dynamical scalar field 6. In 
order to avoid the difficulty with fourth-order derivative terms appeared in f{R) gravity, we 
first transformed the f{R) gravity into the scalar-tensor theory by introducing a scalaron (p. 
This will provide the DCS modified gravity with two scalars, which provides four physically 
propagating degrees of freedom. 

Interestingly, the perturbation equation for the CS scalar 6 is coupled to the odd-parity 
metric perturbation equation, providing a system of two coupled second-order equations, 
while the scalaron (p is coupled to the even-parity perturbation equation. This enables us 
to perform the stability analysis of the Schwarzschild black hole obtained from /(i?)+DCS 
modified gravity theory completely. It was shown that the CS coupling affects the Regge- 
Wheeler equation significantly, while f{R) gravity does not affect the Zerilli equation. It 
turns out that the Schwarzschild black hole is stable against four external perturbations of 
{M, $, Q, ip} if the scalaron is free from the tachyon. 

However, the role of DCS term is limited here because its perturbation SCfj^i, in fl2.16p 
does not involve third-order derivative terms. This higher derivative may appear when 
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the background solution contains a spherically symmetric CS scalar [13]. In this case, 
one could not decouple five massive gravitons successfully because there exists a mixing 
between odd- and even-parity modes, and third-order derivative terms are present. Even 
in the Minkowski background, it is not clear which modes are propagating with their own 
masses. It has been argued that a spacelike vector of v'^ = 8^9 = (0, fl) renders the theory 
free from ghosts and tachyons, while a timelike vector of v'^ = (yU, 0) yields an inconsistent 
quantum theory [23l[2l]. On the contrary, the opposite case is true: the only tachyon- and 
ghost-free model is the one with a timelike vector [25]. Hence, it seems to be a formidable 
task to perform the stability of the Schwarzschild black hole when including the third-order 
derivative terms. 
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Appendix: The explicit forms of twenty perturbation 
equations of E'^ = and Oi = where 



El 



2r5 



2r^k^K{r) + r(-2MV"^ + r^X)Ho{r) + 2i(3M - 2r)r^kHi{r) 
r/(2M2 - r^k^)H2ir) + 2Mr^fK'{r) + r^(5M - 2r)i7^(r) 



3 



+ ^<^ (2M - Mr + fcV^)$(r) + Mr^f^'{r 



Eo = — ^ 
2r2 



2(r - 3M)kf~^K{r) - 2krfH2{r) + Akr'^K'ir) - iXHi{r) 



r / 



2r2 



2A;r 



fc^2^(^) _ 2iMHi{r) + kr^fHiir) - ir^fH[ + ^<l>(r 



^4 = ^-Jy^ 2M(2r - 3M)f-^Ho{r) + 2iMkr^Hi{r) - /(^GM^ - 4Mr 

- r^fX^H2{r) - 2r{QM^ - 7Mr + 2r'^)K'{r) - Mr'^H'^ 
+2ikr^fH[ + rf{QM^ - 7Mr + 2r'^)H'^ - 2r^fK"{r) + rV^^o 



-V"$(r) - J-|(-5M + 2r)$(r) + r(5M - 2r)<l>'(r) + r^-^-''^} 

(r - M)rf"^Hoir) - ikr^Hi{r) - {2M^ - 3Mr + r'^)H2{r) 

+r^fK'{r) - r^H^ + ^| - 2<l>(r) + r<l>'(r) 

6 I 



2r3/ 



2Mrr^Ho{r) - 2ikr^Hi{r) + 2(2M2 + Mr - r^)H2{r) 
+r^{k^r^f-^ + A + 2)A'(r) - 2r2(3M - 2r)ir'(r) - r^H'^ir) 

-r^fH'(r)+r^fK"(r) + 2rf\—r^V"^(r) + i(r^'(r) - $(r^ 



1 

2r 



rf ^Ho{r) - rfHiir) - ^^(r) 
Eg = -Ercotipi, 
Eg = EgsinVi, 
i?io = Ej cos sin ipi, 
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CSC ifi 

2r3 



ir'^fl 2khi{r) + krh[{r) - irh'^ir) \ + (4M - Ar)/io(r) 



Mf\2ip{r) -rip'{r) 



2r^/sin</)i< 2khi{r) + krh[{r) — irh'^{r) > — [rf + 2Mcos2</?i) cscLpihoir) 



O4 



CSC (yCl 

-1 
2^ 



6 r 

+r(csc ipi — X sin (pi)ho{r) — ^—Mf sin < 2?/^(r) — rip^r) 

2tkr'^ho{r) - ikr^h'^ir) + (2r/ + A;V - Xrf)hi{r) - -^kM^{r) 

(f)r 

ikr"^ sin (y9i{2/io(r) — rh^} + \ — rf cos (^1 cot (pi + (r/ + k'^r^) sin (yJi 



62 

+r/(csc</9i — Asiny9i) \hi{r) — -^kM sin (piip{r) 

0r 



CSC cot 
O5 = T7 

Oq = —O5 tan(/Pi, 



ikr^hoir) + r/|2M/ii(r) + r^fh[{r] 



O7 



O5A 



sin (/^i tan Lfi 



Os = O5 sin^ ipi tan <y9i , 
O9 = -O5 sinVi, 
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